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I The expression for the spin susceptibihty x of degenerate quark matter is derived 

o ■ 

' with corrections upto 0{g'^ In g'^). It is shown that at low density, x ^ changes sign 

(N ■ 

^ ■ and turns negative indicating a ferromagnetic phase transition. To this order, we 



o 



also calculate sound velocity ci and incompressibility K with arbitrary spin polar- 



' ization. The estimated values of ci and K show that the equation of state of the 

(N ■ 

polarized matter is stiffer than the unpolarized one. Finally we determine the finite 
Q^. temperature corrections to the exchange energy and derive corresponding results for 

Dh' the spin susceptibility. 
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§ ■ I. INTRODUCTION 

a^ ■ 
o 

One of the active areas of high energy physics research has been exploration of the so 
• 1—1 . 

^ ' called Quantum Chromodynamics (QCD) phase diagram. In particular, with the advent of 

H : 

5r , ultrarelativistic heavy ion beams at RHIC and CERN and with the upcoming facilities of GSI 
where compressed baryonic matter is expected to be produced, such studies have assumed 
special importance. Beside the laboratory experiments, various astrophysical objects like 
neutron stars, quark stars, provide natural sites where many of the theoretical conjectures 
about the various phases of quark matter can be tested. The latter, in the present context, 
is more relevant here, as we study the possibility of para-ferro phase transition in dense 
quark system interacting via one gluon exchange. 

The original idea about para-ferro phase transition in quark matter was proposed recently 
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in [l| where the possibihty of Bloch hke phase transition [2| was studied and it was shown 
that spin polarized quark matter might exist at low density j^. The underlying mechanism 
of such a phase transition is analogous to what was originally proposed for the degenerate 
electron gas Q|. There, for Coulomb interaction, it was shown that the exchange correction 
to the energy is attractive which at the low density wins over the kinetic energy giving rise 
to a ferromagnetic state j^. In [l], a variational calculation has been performed to show 
that it is indeed possible to have a spin polarized quark matter at low density of strange 
quark system, while for light quark it never happens ll]. Similar difference of the light and 
strange quark matter, albeit in a different context, was observed earlier {4]. However, in j^, 
it was shown that both the light and heavy flavor systems can exhibit such phase transitions 
although the critical density for the strange matter is higher than the light quark systems. 
Such investigations, have also been performed in and also in [9, uM where the 

calculation has been extended to include thermal effects. The Bloch like phase transition, 
for strange quark, has also been reconfirmed in 



One shortcoming of all these works including 



has been that the calculations were 



restricted to the Hartree Fock 
termed as correlation energy 



evei and the terms beyond the exchange diagrams, commonly 
1^ . isl . 14, Q, 16] were ignored. Without such corrections. 



however, the calculations are known to remain incomplete as the higher order terms are 
plagued with infrared divergences arising out of the exchange of massless gluons, indicating 
the failure of the naive perturbation series. We know that this problem can be cured by 
reorganizing the perturbation theory where a particular class of diagrams, viz. the bubbles 
are resummed in order to obtain a finite result. Originally, as is well known, this was done 
by Gell-Mann and Brueckner p/7| while calculating the ground state energy of degenerate 



electron gas. The contribution of the bubb 



perturbative nature of the correction 
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es mvo. 
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2l| 



ve terms of 0{g'^ In g^) indicating non- 



In the present work, as announced, we calculate the spin susceptibility (x) of dense 
quark system with corrections due to correlations i.e. containing terms upto O^g^lng"^). 
This requires the knowledge of the ground state energy (GSE) of spin polarized matter with 
inclusion of bubble diagrams. The GSE of the polarized quark matter has been calculated 
only recently in IGj which is the starting point of the present paper. This work is very similar 



to that of Brueckner and Swada 22[ and those of 23 



2J], applied to the case of QCD matter. 



Unlike, degenerate electron gas, however, we have both the electric and magnetic interactions 



3 



and the calculation is performed relativistically, while the non-relativistic results appear as 
a limit. 

The spin susceptibility Xi for quark matter upto 0{g'^) has already been calculated in 
Ref. |l| which we only briefly discuss. Subsequently, the non-fermi liquid corrections to x 
has also been studied in [oi These studies provide further motivation to undertake 

the present endeavor to include correlation corrections, without which, as mentioned al- 
ready, the perturbative evaluation of x remains incomplete. In addition, we also calculate 
incompressibility and sound velocity for spin polarized quark matter with corrections due to 
correlations which involve evaluation of single particle energy at the Fermi surface. These 
quantities are of special interests for applications to astrophysics. Moreover, we also evalu- 
ate the exchange energy density at non-zero temperature and determine the corresponding 
corrections to the spin susceptibility. 

The plan of the paper is as follows. In Sec. II we calculate spin susceptibility with 
correlation correction for degenerate quark matter. Analytic expressions are presented both 
in ultra-relativistic (UR) and non-relativistic (NR) limit. In Sec. Ill, we evaluate exchange 
energy density and spin-susceptibility at non-zero temperature. In Sec. IV we summarize 
and conclude. Detailed expressions of the intermediate expressions, from which x is derived, 
have been relegated to the Appendix. 



II. SPIN SUSCEPTIBILITY 



The spin susceptibility of quark matter is determined by the change in energy of the 
system as its spins are polarized 22]. We introduce a polarization parameter ^ = (n^ — 
n~)/nq with the condition < < 1, where n'^ and n~ correspond to densities of spin-up 
and spin-down quarks respectively, and Uq = +n~ denotes total quark density. The Fermi 
momenta in the spin-polarized quark matter then are = P/(l+0^^^ Pj = 
where pj = (vr^rag)^/^, is the Fermi momentum of the unpolarized matter = 0). In the 
small ^ limit, the ground state energy behaves like jl| 



(1) 
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Here, Ps 



defined to be the spin stiffness constant in analogy with ly, |2l|]. The 



spin susceptibihty x is proportional to the inverse of the spin stiffness, mathematically 
X = 2/3~^ [25|. It is to be noted that in Eq.([l]), the first term corresponds to unpolarized 
matter energy. 

Now, the leading contributions to the ground state energy are given by the three terms 



viz. kinetic, exchange and correlation energy density 



16| i.e. 



^kin ~l~ ^ex ~l~ ^corr • 



The total kinetic energy density for spin-up and spin-down quark becomes 



m 



(2) 



kin 



167r2 



-nig In 



1 + 0'/Vp?(1 + 0^/=^ + [2pJ(l + 0'/' + ml] 



m„ 



rrin 



(3) 



where rriq is the quark mass. 

The exchange energy density E^x have been calculated in ref. within Fermi liquid 
theory approach. One can also directly evaluate the two loop diagram [l| to obtain 



eL 
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(4) 
(5) 



where f^^, and , stands for non-flip (s = s') and flip (s 



—s') forward scattering 

amplitude given in [ill, ll6| . Here, E'ex = E^^+E^^ can be estimated numerically. However, 
analytical evaluation of these integrals is possible in the ultra-relativistic and non-relativistic 
limits as reported in [l, [ill, Q. 

The next higher order correction to t he g round state energy beyond the exchange term is 



the correlation energy Ecorr [121, Q, Q, ll5|]. The detailed calculation of correlation energy 
for spin polarized matter have been derived in 1^ which we quote here: 



E ~ 

-'-'corr — 



1 1 



(27r)3 2 



n/2 



sin^ OEdOE { H 



In 



H, 



2H| 



In 



H. 



(6) 



with 9e = tan ^{\k\/kQ). The relevant Hl and are determined to be [l 
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III 



III 



3_ 

47r2 

87r2 



E- 



sin^ 



cot 6^^ -1 / Q /, \ 
1 ^ — tan \v^ta,iiUE) 



cot 



s=± 



cot 6'r 



The spin susceptibihty is given by 



cot^ Be 



tan ^ {y'j tan^g) 



(7) 

(8) 



1 d^E{0 



We have = Xkln + XfTx + XiJorr- The kinetic and exchange contribution have been 
evaluated in ref . [l| , is given by 



(10) 



6p2 



3P/ 



j>/£:/ — In 



fpf + ^f\ 


r 







1 + 



2m„ 



Hpf 



To determine the correlation correction to spin susceptibility, we expand curly braces 
terms of Eq.([2]) in powers of the polarization parameter ^, which gives 



In 



-2 
''f 



2m 



In 



{AoL + Bot) + e{AiL + Bit) + 



(12) 



Here, Aql and i3oT correspond to unpolarized matter term and the detailed expressions of 
Ail and Bit are given in the Appendix. Xcorr 



Xcorr 



corr 

2 9^ 



1 1 



(27r)3 2,o 



tt/2 



sin^ 9Ed9E{AiL + B,T) 



(13) 



From the above expression and with the help of the expression presented in the Appendix, 
Xcorr can be estimated numerically. Results for the two limiting cases however can be 
obtained analytically as we present in the following two sub-sections. 
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A. Ultra-relativistic limit 



In the ultra-relativistic limit, the kinetic, exchange and correlation energies are 



rpur 



3pj 



[{I + if' + {I - 0'"] 



J^// [(1 + er/^ + (1 - er/^ + 2(1 - en , 



With the help of Eq.([T]), each energy contribution to the susceptibility is 



(14) 



Xe 



A. 

67r^ 



\corr 



'367r4 

4 4 

g Pf 



(Inr, -0.286). 



23047r6-"'^ ^^^^ 

with r, = g^y/'"^. From Eq.(fT5D. sum of all the contribution to the susceptibility can 
be written as 16 1 



X" 



xp[i- 



g' 



g^ 



67r2 3847r4 
where xp is the non-interacting susceptibility 



(Inr, - 0.286)]-S 



(16) 



23 



24|. 



B. Non-relativistic limit 

Now, we go to the non-relativistic limit to calculat e sp in-susceptibility in order to compare 
our results with those of dense electron gas 



22 



23 



24 



26| interacting via, static Coulomb 
potential. In this limit, kinetic and exchange energy densities are jl|, [ill, 16|: 



rpnr 



3p 



f 



kin 



207r2m, 



[(1 + 0^/^ + (1-0^^1 



g^ 



4/31 



(17) 
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The contribution to the susceptibihty from kinetic and exchange energy density yields 



Xe 



(18) 



We want to calculate the contribution to the spin-susceptibility beyond the exchange 
correction. For this we first evaluate the correlation energy in this limit. 
The dominant contribution to the correlation energy is found to be, 



n^mg Jxi/2 k 



— / xdx^/(s) 



s=± 



l--lnf^)]E^(l--')' (19) 

^ s'=± 



where A = {g'^rrig) / {8np j) , f{s = ±) = (1 ±^)^/'^, x = x'^f{s), x^ = (konig) / (pjk) and 
k' = k/pf. For s = s' one obtains: 



g^lng'^l 3 
-mgp^(l -ln2), 



(20) 



(27r)6 3 

Note that, here, the correlation energy is independent of spin-polarization ^. It is seen, 
that for the spin parallel interactions .^-dependent terms contribute with opposite sign and 
cancels each other. For s = —s', the integral on x takes the form 



xdx < (1 



3^) 



l-^a:(l-^01n 



x(l -10 + 1 



x(l - iO - 1 



^[i-x(i + -o] + (e--o 



Expanding In in terms of C, and retain upto Ol^"^) we have 



f211 



Using Eq.([in]), dni) and 



3 

we have 



:i-ln2) 



6^ 



(22) 



- 1287r6 3 ' 



:i-ln2) 



6^ 



(23) 
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It is to be mentioned that similar expressions for degenerate electron gas interacting via. 
static Coulomb potential can be found in ref. 26|. From Eg. (1201) and fl23p . it is clear that 
spin anti-parallel states are attractive in contrast to the parallel states due to Pauli exclusion 
principle. In this limit the correlation contribution to the susceptibility is found to be, 



X.corr 



9 In 3 



(24) 



The total susceptibility is given by 



r 



Xp 



1 - 



3847r4 pj 



(25) 




In Fig.([T]) we plot inverse spin susceptibility which is valid for all the kinematic regimes. 
It shows changes its sign at the density ~ 0.12fm~'^ without correlation correction and 
when we include the correlation effect its sign changes at ~ O.lfm"^. This is equivalent to 
what happens to the ground state energy as a function of ^. It is needless to mention that 
this change of sign correspond to the para-ferrophase transition in dense quark system. The 



parameter set used here are same as those of 
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C. Incompressibility and Sound velocity 

Once we have the expressions for the total energy density, the incompressibihty (K) and 
sound velocity (ci) can be determined. The incompressibility K is defined by the second 
derivative of the total energy density with respect to the number density Uq, which is given 



by 



ill 



K ^ 9n,g. (26) 
Since, there are two Fermi surfaces corresponding to spin-up (+) and spin-down (— ) 



states, such that E = E(n^, ). We have 



ill 



dE _ dE <9n+ dE dn- 



dug dw^ duq dn^ dug 

= ^[(1+0/^^ + (1-0/^1- (27) 

The single particle energy at the Fermi surface or chemical potential of spin-up quark 
turns out to be 



9' ( .Pf 



7687r4 



Pf_ 



(28) 



Similarly, /i"'"'' can be obtained by replacing with pj in Eg. (1251) . In ref. ll|, chemical 
potential is determined within Fermi liquid theory approach upto 0{g'^). We, however, here 
calculate jji^ with different approach upto 0{g'^\n.g'^). 

Using Eq. fl27|) and Eq. fl28l) . the incompressibihty becomes 



= \pf {[(1 + + (1 - if'] + Yf^[(i + if" + (1 - if" + 2(1 - er'] 

e?"]]. (29) 



+?i?^[(i + 0'/' + (i-0'/' 



2 1 



7687r4 J 

Another interesting quantity would be to calculate the first sound velocity which is given 
by the first derivative of pressure with respect to energy density. Mathematically 
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(1 



^)^<? dn- 



;i + 0/"+ + (1 - 



(30) 



From Eq. (l28|) . we have 



3p2(l + ^)2/3 



9' 



7687r4 



127r2 

^)2/3_(. 



;i + ^)2/3_(l_^)2/3 



;i + e)2/3 







2/3 



^1+^)2/3 f (31) 

The second and last term in the curly braces corresponds to exchange and correlation 
contribution respectively. Similarly, dfi~/dn~ can be obtained by replacing ^ with — ^. 
Using n^, js^ and dfi^/dn"^, we calculate the sound velocity in terms of ^. Numerically, for 



unpolarized matter ci = 0.46, while for complete polarized matter ci 
the causal value l/\/3 = 0.57 at the high density limit. 



0.54, which is below 




FIG. 2: Incompressibility K in the spin polarized quark matter. 



In Fig ([2]), we plot the density dependencies of the incompressibility with correlation 
correction. This shows for higher value of the order parameter C,, the incompressibility 
becomes higher for the same value of density. Thus numerical values of incompressibility 
and sound velocity shows that equation of state for polarized quark matter is stiffer than 
the unpolarized one jll |. 
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III. SUSCEPTIBILITY AT NON-ZERO TEMPERATURE 



In this section we calculate the exchange energy density E^.^ at low-temperature (T << 
for which we replace 0{p^ — of Eqs.(j4ll5]) with proper Fermi distribution function. 
In the ultra-relativistic limit, the angular averaged interaction parameter is given by 111 ] 



£ur 
J pp' 



9' 



dl^i f dn 



1 + (p ■ s){p' ■ s) 



9pp' J Att J An 
The spin non-flip contribution to the exchange energy density is 



(32) 




dp dp nf s/TTi\ s /rr\ 



(27r)3 (27r)3- 



-E 



487r2 



2„2 



2/3 



:i-o 



2/3 



(33) 



Here ^p[pi){T) is the Fermi distribution function. 

Similarly, El^ can be evaluated. The total -E"^' at low temperature is found to be 



ex — 



9' 4 



327r 



247r2 



l+^)4/3+(l_^)4/3^2(l-e')'/' 

1 + 0'/' + (1 - 0'/' 



(34) 



The kinetic energy density can be written as 



TP' 



?,p) 



kin 



1 + 0'/' + (1 - 



4/3 



+ 



3r 



2„2 
Pf 



[(1 + 0^/^ + (1 - 0^/^ • 



2/31 



(35) 



From Eq.([T]) each energy contribution to the susceptibility is 



Xe 



P} ^ 



2rp2 



P) 



9Yf , TT^T^ 



367r4 



?>p) 



(36) 



It is to be noted that the T independent terms of the above expressions are identical with 
those Eqs. fll4tiT^ . Thus the susceptibility at non-zero temperature is given by 
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Xp 



67r2 V" ' 3j9^ 

In the non-relativistic limit the interaction parameter takes the following form 



(37) 



rnr 
J pp' 



9 



1 + s ■ s' 



(38) 



For spin anti-parallel interaction s = —s', then f^^, = 0. Thus the contribution due to the 
scattering of quarks with unlike spin states vanishes and the dominant contribution to energy 
density comes from the parallel spin states (s = s'). Performing the angular integration of 
Eq.(jl]), the exchange energy density upto term (9(T^) becomes 



rrinr 



47r4 



Y^jpdV n;{T) jp'dp' n;,{T)\n 

^)4/3 + (1 _ ^)4/3 



s=± 



p + p' 



p — pi 



The kinetic energy density is found to be 



l-^T^m,pf[{l + iY/^ + {l-if/^ 



3p^ 



f 



207r2m„ L 



1 + 0'/' + (1-0 



5/3 



+ 



1 + ^)2/3 + (1 _ ^) 



2/3 



(39) 



(40) 



Separate contribution from kinetic and exchange energy to susceptibility becomes 



Xkin 



Xe 



p) 



1 - 



3pj 



IM*V 2pj 
Thus, at low temperature the susceptibihty turns out to be 



(41) 



X" 



Xp 



S-K^Pf 



6p^f 



(42) 
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IV. SUMMARY AND CONCLUSION 



In this work we have derived the spin susceptibihty for degenerate quark matter with 
corrections due to correlation contributions. Analytic expressions for susceptibility are also 
derived both in the ultra- relativistic and the non-relativistic limit. It is observed that at 
low density susceptibility changes sign and becomes negative suggesting the possibility of 
ferromagnetic phase transition. In addition, we also derive single particle energy, sound 
velocity and incompressibility upto 0{g^ In (7^). As far as the equation of state is concerned, 
in the present model, we find that equation of state for polarized matter is stiffer than that 
of unpolarized one. We also determine the exchange energy and susceptibility at non-zero 
temperature of the spin polarized quark matter. 
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V. APPENDIX 



To calculate the correlation contribution to the spin susceptibility we have from Eq. flT^ 



g^pj- sec^ 9e csc^ 9e 

'^^^ ^ ~ 11527r%3(r7i2+p2sec2 OeY 

.2 



X In 



' g"^ csc^ 6e 

2 



Pf — Ef cotOE tan ^(vftanOE 
tan 9 e) {p^f + m^cos^ 9 E 



X |64£:^cos^6'£;tan ^{vf 

- 2pfej sin 29e tan~^{vf tan 9E){l2ml + SlmJpJ 



4^2 



m^g{Aml + 5pJ) cos46'£; + 68m^pJ + Aml{Am^g + Wm^pj + 7pj) cos26'£; 

' 29mJpJ + mJ(2mJ + 3p5)cos4^i5 + 36mJp} 



■32^6} 



+ 4pj5/ sin^ ^£;[6mg + 2ymgPj , , 
+ 4mJ(2mJ + 4mJpJ + 3p)) cos 29 e + 16|>5] | 



(yf^Pj cot^ 9e CSC* 9e 



_ y ^/ ^ ^ 
" 11527r%3(m2 cos2 ^^j+p^) 

^2 cot CSC2 r _i 2 2/1 , 2n • o/i 

^ ^ 2 tan (17/ tan6'£;j(mg cos 9e + Pj) — Pf£fSin29E 



X In 



X 



I - 324 tan"^(t;/ tan 6'£;)(mg cos^ 9e + p}f 

QJ2^ r™4 , ^4 , ^JpJ(i + cos2^B)]sin2 2^i5 + 2pjtan-^(t;jtan^E)sin2^£; 



- ^p]ef[m\+p) + mlp]^ 
X [SmJ + 31m^pJ + m^pj cos A9e + 36m^p) + 4mJ(2mJ + 4m^pJ + 3pJ) cos 29 e + 



(43) 



16//]}. 
(44) 



with Vf = pf /Ef. 
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